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Abstract. We present the results of a series of optical, UV, X-ray and 7-ray observations of the BL Lac object 
S50716+714 carried out by the Swift and AGILE satellites in late 2007 when this blazar was flaring close to its 
historical maximum at optical frequencies. We have found that the optical through soft X-ray emission, likely 
due to Synchrotron radiation, was highly variable and displayed a different behavior in the optical UV and soft 
X-ray bands. The 4-10 keV flux, most probably dominated by the inverse Compton component, remained instead 
constant. The counting statistics in the relatively short AGILE GRID observation was low and consistent with a 
constant 7-ray flux at a level similar to the maximum observed by EGRET. An estimate of the 7-ray spectral slope 
gives a value of the photon index that is close to 2 suggesting that the peak of the inverse Compton component 
in the Spectral Energy Distribution (SED) is within the AGILE energy band. The different variability behavior 
observed in different parts of the SED exclude interpretations predicting highly correlated flux variability like 
changes of the beaming factor or of the magnetic field in simple SSC scenarios. The observed SED changes may 
instead be interpreted as due to the sum of two SSC components, one of which is constant while the other is 
variable and with a systematically higher synchrotron peak energy. 



1. Introduction 

Blazars are a rare class of radio loud AGN known to 
display extreme properties like irregular, often large and 
rapid variability, apparent super-luminal motion, flat ra- 
dio spectrum, large and variable polarization. Blazars are 
thought to be objects emitting radiation over the entire 
electromagnetic spectrum from a relativistic jet that is 
viewed closely along the line o f sight thus causing st rong 
relativistic amplification (e.g., lUrrv fc Padovanil (|l995l )). 
Blazars come in two flavors: Flat Spectrum Radio Quasars 
(FSRQs, or blazars of the Q type in the nomenclature of 



Giommi et al 



2008) , and BL Lac objects, (or blazars of 



Giommi et al. 2008) depending on whether 



the B type in 

their optical spectrum does or does not show the strong 
and broad emission lines that are typical of QSOs. The 
broad-band continuum emission of both types of blazars is 
normally associated to synchrotron and inverse Compton 
emission, which is expected to produce the two-bump 
shaped Spectral Energy Distribution (SED) that is always 
observed in these sources. The synchrotron component 
usually peaks [in a Log{vf{v)) — Log(v) representation] 



in the infra-red or in the optical band but is some cases, 
almost exclusively in Blazars of the BL Lac type, the syn- 
chrotron emission can peak at much higher energies, up 
to the hard X-ray band. 

S5 0716+714 is an optically bright and well stud- 
ied BL Lac known to be highly variable at all frequen- 
cies (e.g. Wagner et al. 19961 : iGiommi et aL 1999f) and on 
time scales ranging from years to a few th o usand seconds 
(e.g. iRaiteri et aL 2003 : Nesci et "all l2005t iFoschini et al. 
2006). Within the classification of IPadovani fc Giommi 
(|l995tl S5 0716+714 is an intermediate object (IBL) be- 
tween Low Energy Peaked (LBL) and High Energy Peaked 
(HBL) blazars since its synchrotron component reaches 
the soft X-ray and the Inverse Com pton component has 



been detected in the hard- X-rays ( Giommi et al. 19991 : 



Tagliaferri et al.ll2003l ). Despite the many attempts to de- 
termine its redshift no positive detection of any feature 
in its optical spectrum has ever been reported. This, com- 
bined with the lack of any host galaxy sets a lower limit of 
0.52 to its redshift (jSbarufatti et alj|2005h . S5 0716+714 
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was detected several times by the EGRET 7-ray de tector 
of the Compton Observatory ( Hartman et al. 1999f ). 

In this Letter we present a series of AGILE and Swift 
ToO observations of S5 0716+714 carried out between 23 
October and 13 November 2007 following an AGILE de- 
tection of this source with a large 7-ray flux in September 
2007 (Giuliani et al. 2007, ATEL 1221) and an alert issued 
in late October 2007 by an optical monitoring campaign 
reporting that this source w as undergoing a large optical 
flare (see IVillata et al.ll2008l for details). 

2. AGILE observations and data analysis 

The AGILE observatory (jTavani et al . 2008) is a new tech- 
nology satellite of the Italian Space Agency (ASI) dedi- 
cated to 7-ray and X-ray astrophysics. AGILE was suc- 
cessfully launched on 23 April 2007 into a low Earth, 
nearly equatorial, orbit which ensures a low particle back- 
ground and short passes through the South Atlantic 
Anomaly. AGILE is capable of observing cosmic sources 
simultaneously in the 18 - 60 kcV and 30 M eV - 50 GeV 
energ y bands with the Super AGILE (SA), ( Feroci et alj 
20071 ) and the Gamma Ray Imaging Detector (GRID) 
respecti vely. The GRID is a combination of a Silicon 
Tracker dPrest et al.ll2003t iBarbiellini et al.ll200lh . a Mini- 
Calorim eter dLabanti et al.|2006l) and an Anti-coincidence 
System (jPerotti et alJ^OOaT " 

On 23 October 2007, S5 0716+714 was reported by 
optical observer s to be flaring at a flux close to its histor- 
ical maximum ([Villata et all 120081 ) . Despite AGILE was 



Table 1. Results of AGILE-GRID analysis 



still undergoing its in-flight calibration program, on 24 
October the satellite began a ToO maneuver to bring 
S5 0716+714 in its very large field of view (~ 60 deg ra- 
dius). The ToO lasted until 1 November 2007 11:50 UTC, 
for a total pointing duration of ~ 8 days. 

All GRID data are normally processed using the sci- 
entific data reduction software developed by the AGILE 
instrument teams and integrated into an automatic 
pipeline system developed at the ASI Science Data Center 
(ASDC). The 8 days of S5 0716+714 October ToO data 
have been processed with the most recent software and 
in-flight calibrations available at the time of writing 
(BUILD-GRID 14). The AGILE software was run creat- 
ing counts, exposure, and galactic background 7-ray maps 
with a bin-size of 0.3 x 0.3 deg. To reduce the parti- 
cle background contamination only events tagged as con- 
firmed 7-ray events were selected (F4 filtercode=5), All 
the 7-ray events whose reconstructed directions form an- 
gles with the satellite-Earth vector smaller than 80 deg 
(albrad=80) were also rejected, thus reducing the 7-ray 
Earth Albedo contamination. We created maps in the full 
AGILE energy band E > 100 MeV, and in the two sepa- 
rate energy bands: 100-400 MeV and 400-1000 MeV. The 
AGILE Maximum Likelihood procedure gives the average 
flux values reported in Table Q] 

A preliminary spectral fitting based on a least square 
fit to the flux estimated in different bands gives a photon 
spectral index of ~ 1.9 + 0.5. The low counting statistics in 
the rather short (8 days) AGILE-GRID observation does 



Energy band 


Count rate 


Net counts 


Vts 




(ph cm" 2 s _1 ) 






> 100 MeV 


(4.9+1.1) x 10"' 


48+11 


6 


100-400 MeV 


(2.3+0.9) xl0~ 7 


28+10 


3 


400-1000 MeV 


(6+2) xlO" 8 


11+4 


4 



not allow us to detect small amplitude variability. The 
overall flux is consistent with a constant value or with 
variability of magnitude less than a factor of «3. A de- 
scription of the S5 0716+714 GRID observati ons and fur- 
ther d etails of the data analysis are reported in lChen et al.l 
(|2008h . S5 0716+714 is too faint to be detected by the 
Super- AGILE X-ray detector. 



3. Swift observations 

Swift performed 10 short ToO observations of 
S5 0716+714 between 23 October and 13 November 
2007 as detailed in Table [2] The observations were carried 
out using all three Sw ift on-board experim ents: the 



X-ray Telescope (XRT, iBurrows et ail |2005|) the UV 



and Optical Telescope (UVOT, Roming et al 



2005h an d 



the Burst Alert Telescope (BAT. iBarthelmv et al.1 120051 ) 



However, the hard-X-ray flux of S5 0716+714 is below 
the sensitivity of the BAT instrument for short exposure 
and therefore the data from this instrument will not be 
used here. 

Table 2. Log of Swift observations in October-November 
2007 



Observation date 



XRT exposure 
seconds 



Count rate 
ct/s 0.3-10 keV 



23-Oct-2007 


2443 


0.47+0.015 


24-Oct-2007 


1991 


0.29+0.013 


25-Oct-2007 


2698 


0.37+0.013 


26-Oct-2007 


2199 


0.34+0.013 


27-Oct-2007 


1961 


0.72+0.030 


28-Oct-2007 


1758 


0.33+0.015 


03-Nov-2007 


2801 


0.90+0.030 


06-Nov-2007 


1923 


0.28+0.013 


09-Nov-2007 


2003 


0.39+0.015 


13-Nov-2007 


2382 


0.27+0.012 



3.1. UVOT data analysis 

During all observations UVOT produced a series of images 
in each of the lenticular filters (V, B, U, UVW1, UVM2, 
and UVW2). Photometry of S5 0716+714 was performed 
using the standard UVOT software distributed within the 
HEAsoft 6.3.2 package and the calibration included in the 
latest release of the "Calibration Database" . Counts were 
extracted from aperture of 5" radius for all filters and 
converted to fluxes using the standard zero points. The 
fluxes were then de-reddened us ing a value for E(B — V) of 
0.031 mag ([Schlegel et al.lll998h with A X /E(B - V) ratios 
calculated for UVO T filters using the m ean interstellar ex- 
tinction curve from Fitzpatrick ( 19991 ) . No variability was 
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detected within single exposures in any filte r. Part of the 



UVOT data presented here also appear in iVillata et al 



(2008) who reported results fully consistent with ours. 



3.2. XRT data analysis 

All XRT observations were carried using the instrument in 
Photon Counting readout mode, which provides maximum 
sensitivity but is affected by photon pile up for count- 
rates larger than w 0.5 cts/s. The data were reduced using 
the XRTDAS software developed at the ASI/ASDC and 
distributed within the NASA/HE AS ARC HEAsoft 6.3.2 
package. Photons were selected with grades in the range 0- 
12 and default screening parameters were used to produce 
level 2 cleaned event files. Spectral data were extracted 
in a circular region with a radius of 20 pixels. To avoid 
pile-up problems, when the source count rate was larger 
than 0.5 cts/s, the spectral data were extracted in an an- 
nular region with inner and outer radii of 3 and 20 pixels 
respectively The background was estimated in a nearby 
source-free circular region of 50 pixels radius. In order to 
use x 2 statistics, spectra were rebinned to include at least 
20 photons in each energy channel. For our spectral fitting 
we used the XSPEC 11.3 spectral analysis package. Given 
the very short exposures (typically w 2,000 seconds, see 
Tab [2]) the limited photon statistics allowed us to fit the 
data to a simple power law spectral model with low energy 
absorption (Nh) fixed to to the Galactic va l ue of 3.8 x 10 20 
cm~ 2 as estimated by Dickey fc Lockman ( 1990h . 
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Fig. 1. The light curves of S50716+714 in the soft XRT 
band (0.3-4 keV, top panel), in the hard XRT band (4-10 
keV, second panel from top), in the UVOT U filter (third 
panel) and in the AGILE-GRID detector (bottom panel). 
The vertical axis of the Swift light curves is logarithmic 
with one decade range for all cases, showing that varia- 
tions are largest in the soft X-ray band. 



The results are reported in Tab. [3] where column 1 
gives the observation date, column 2 gives the flux in the 
2-10 keV band, column 3 gives the power law photon index 
[fiph cm~ 2 s~ 1 ) oc E~ r ] and ler error, and column 4 gives 



the reduced % 2 and the number of degrees of freedom. In 
all cases the average photon index is steeper than 2, indi- 
cating that the 0.3-10 keV X-ray spectrum is dominated 
by synchrotron radiation. However, from the X-ray light 
curves of Fig. 1 clearly show that the soft (0.3-4 keV) X- 
ray flux varies much more than the hard (4-10 keV) X-ray 
flux, suggesting that the high energy part of the spectrum 
is mostly due the rise of the inverse Compton component 
as is also apparent in the broad-band SED of Fig. 2. 

Table 3. Results of XSPEC fits. Power law model with 
Ah fixed to the Galactic value of 3.8 x 10 20 cm" 2 . 



Observation 


2-10 keV flux 


Spectral slope 


Reduced 


date 


erg cm 2 


s" 1 


r 


X 2 /d.o.f. 


23-Oct-2007 


5.9x10" 


1'2 


2.3 ±0.05 


0.74/40 


24-Oct-2007 


4.5x10" 


12 


2.2 ±0.08 


0.55/18 


25-Oct-2007 


4.8x10" 


12 


2.2 ±0.05 


0.81/35 


26-Oct-2007 


5.1x10" 


12 


2.1 ±0.07 


1.37/26 


27-Oct-2007 


7.6x10" 


12 


2.4 ±0.10 


0.68/32 


28-Oct-2007 


5.6x10" 


12 


2.2 ±0.07 


0.87/19 


03-Nov-2007 


8.0x10" 


12 


2.6 ±0.07 


0.98/54 


06-Nov-2007 


3.9x10" 


12 


2.4 ±0.08 


1.28/17 


09-Nov-2007 


5.6x10" 


12 


2.4 ±0.07 


1.08/24 


13-Nov-2007 


3.7x10" 


12 


2.4 ±0.07 


0.92/21 



4. Discussion 

We have presented a series of simultaneous ToO observa- 
tions of the blazar S5 0716±714 carried out by the Swift 
and AGILE satellites following the detection of this so urce 
in a state of large flaring activity (jVillata et al. l l2008lL 
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Fig. 2. The SED of S50716±714 built combining our 
data with s i multan eous optical and radio data from 
IVillata et al. ( 2008f ) (color) and with archival multi- 
frequency data (light grey). The dashed and dotted lines 
represent the SSC components described in the text. 



The light-curves shown in Fig. 1 show strong variabil- 
ity (up to a factor rj 4) in the soft-X-ray flux, more mod- 
erate variations (less than a factor 2) at optical/UV fre- 
quencies, and an approximately constant hard X-ray (5-10 
keV) flux (see also Fig. 2). A similar behavior was seen 
during a BeppoSAX observation of S5 0716±714 carried 
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out in October 2000 ( Tagliaferri et al.ll2003t Foschini et al. ken power law spectrum for both particle components: 



2006; iFerrero et alj 120061 ). The low counting statistics ac 
cumulated during the relatively short (~ 8 days) AGILE 
observation does not allow to detect small amplitude vari- 
ability. The count-rate is consistent with a constant flux 
or with variations lower than a factor of « 3. 

The Spectral Energy Distribution of S5 0716+714 built 
with all the Swift and AGILE data is shown in Fig. 2 where 
we see that our flux measurements (colored filled circles) 
are much larger than those from archival data (light grey 
points). Large intensity and spectral changes close to the 
peak and at the end of the synchrotron component, are 
clearly present while the rise of the inverse Compton com- 
ponent above approximately 4 KeV remained constant 
over the entire observing period. The AGILE GRID data, 
integrated between 23 October and 1 November, is shown 
as three points between wlOO MeV and 1 GeV which have 
been derived from the fluxes reported in table [TJ The 
intensity is close to the maximum observed by EGRET 
(grey points). This flaring behavior is quite different from 
that seen in 3C454.3 in April/May 2005 where the in- 
verse Compton componen t varied by a large factor in a few 
days ( Giommi et al J 2006h . This demonstrates that flaring 
events in blazars may be widely different and each should 
ideally be followed with detailed observations covering the 
dynamical times scales of all the emission components. 
The observed spectral behavior of S5 0716+714 cannot be 
the result of simple changes in the magnetic field or of the 
beaming factor in a one-zone homogeneous SSC model. In 
this case the emission in different bands is expected to vary 
in a highly correlated way, contrary to what is observed in 
the optical and X-ray data of Fig. 1 and in the we ll sam- 
pled simultaneous light curve of lVillata et alJ (|2008l ) where 
the R-band variations are considerably smaller then in the 
soft X-ray band. This behavior would also challenge exter- 
nal C ompton model (jSikora et al ] |l994t iTagliaferri et all 
2003f ) in which variable high energy emission should result 
from a single electron population accelerated by distur- 
bances in the jet. 

The different flux variations observed at optical and 
X-ray frequencies could instead be explained by the pres- 
ence of two SSC components, one of which is constant over 
the entire observing period, while the second one is highly 
variable and possibly due to a secondary blob of relativis- 
tic plasma including fresh and more energetic particles 
(see Fig. 2). However, even this scenario cannot repro- 
duce exactly the light curves and spectral changes in all 
energy bands as this would require a very good knowledge 
of all the acceleration events and a proper treatment of the 
energy-dependent loss mechanisms. This can be achieved 
only through a continuous monitoring of the source while 
our optical, UV and X-ray observations are just snapshots 
of w 2,000 seconds taken approximately at one day inter- 
vals. Fig. 2 shows the two SSC components as dashed and 
dotted lines, while the continuous line is their sum. For 
their evaluation we have assumed a SSC model with bro- 



n e (i) = fc 7 Pl 7 < 76 

n e {l) = k lb Pl {l /lb)~ P2 7>76i 



(1) 
(2) 



with constant density within a spherical region of radius 
r and with 7 in the range 10 2 — 10 6 . The parameters of 
the "stationary" component, which reproduces the "min- 
imum" fluxes at UV and X-ray frequencies, are p\ = 1.7 
and P2 = 4.2 at low and high energies respectively, and 
the break Lorentz factor is % = 10 3 25 ; for the second 
component similar parameters are assumed (pi = 1.1 
and pi — 4.05), with a larger break Lorentz factor of 
7b = 10 3 ' 6 , which produces a stronger emission at high 
energies. The magnetic field and the beaming factor for 
the two components are assumed to be the same (B = 1.1 
G and 5 = 20), while the first component has an higher 
normalization and a larger radius (fco = 631. cm~ 3 and 
r = 0.010 pc) w.r.t. the second component (kg = 8.91 
cm~ 3 and r = 0.005 pc). B and r values are co-moving 
with the emitting regions. 

Finally we note that during a recent large flare 
dGiommi et al.ll2008l) S 5 0716+714 was detected at TeV 
energies ( Teshima 2008[) despite its likely hi gh redshift. If 
the limit of z > 0.52 (|Sbarufatti et al.|[2005h is confirmed 
this detection would imply that the Universe is much more 
transparent to TeV photons than previously claimed. 
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